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Abstract 

We study the impacts of reheating temperature on the inflationary predictions of the 
spectral index and tensor-to-scalar ratio. Assuming sinusoidal oscillations and that re¬ 
heating process is very fast, the reheating temperature can be constrained for sinusoidal 
oscillation within a factor of 10 - 100 or even better with the prospect of future observa¬ 
tions. Beyond this, we find that the predictions can also be insensitive to the reheating 
temperature in certain models, including Higgs inflation. 



1 Introduction 


Recent observational advances place more stringent bounds on the predictions of inflation, for 
example the Planck mission [T] . When the modes relevant for the cosmic microwave background 
(CMB) observations exit the horizon during inflation, the slow-roll approximation is supposed 
to be working very well. Consequently, the properties of the primordial curvature perturbation, 
which after inflation becomes the seed of the observed temperature fluctuations in the CMB as 
well as large scale distribution of galaxies, are determined by the form of the inflaton potential in 
simple canonical models [2]. For example, the spectral index n-ji and the tensor-to-scalar ratio r, 
see ([^ and ([^ respectively. More precise estimates of them, however, require full understanding 
of the thermal history of the universe beyond the end of the inflationary phase. The largest 
terra incognita is the stage between the end of inflation and the beginning of the standard 
radiation dominated epoch, during which the energy in the inflaton sector is transferred to the 
standard model sector [3]. The detail of this reheating process is much less understood, as the 
huge range the reheating temperature Trh may lie within indicates, 1 MeV < Trh < 10 Ge\Q 

Nevertheless, like single field slow-roll inflation, there is a simple, minimal scenario of reheat¬ 
ing 0 . In this scenario, the inflaton oscillates around the minimum of its effective potential 
after inflation ends, with a decreasing amplitude due to the damping by the cosmic expan¬ 
sion, until the expansion rate H falls down below the decay rate T of the inflaton. Then the 
energy of the inflaton sector is transferred into the standard model sector efficiently and the 
universe becomes dominated by relativistic particles, corresponding to a reheating temperature 
Trh ~ A/mpiT 0. Therefore in this minimal scenario of reheating, there are two parameters 
needed for deriving the duration of reheating, an effective equation of state Wrh to describe 
the expansion rate during the oscillating phase of the inflaton, and the decay rate T to deter¬ 
mine when reheating happens. The constraint to the reheating temperature can be derived 
from the observational constraints of n-ji and r once w^h is known. Given the potential form, 
Wrh can be derived, for instance as in the case of power-law chaotic inflation [7j. For a typical 
case where the local minimum is quadratic and tCrh ~ 0, in a universal class of model where uti 
is inversely proportional to the number of e-folds during inflation, we can constrain Frh within 
a factor of 100 or so. With future precision experiments, this constraint can be further reduced 
to a factor of 0{1) or even better. 

If the effective potential is not quadratic at the minimum where Wrh 7^ 0, the above esti¬ 
mation for the reheating temperature may not be valid. In an extreme case where tCrh = 1/3, 
the observational data cannot give any constraint on T^h. Equivalently, it is secure to say that 
the predictions for n^i and r become robust and independent of the detail of reheating in this 
case if we assume the minimal scenario of reheating. We find that the 0^ model, even with a 
non-minimal coupling to gravity such as the case of Higgs inflation, has this property. 

This article is organized as follows. In Section we revisit how the post-inflationary 
evolution, particularly reheating, leads to a theoretical uncertainty in the number of e-folds of 
expansion during inflation, AW- This theoretical uncertainty is then translated into those of 
the inflationary model predictions such as n-ji and r, which depend on the effective equation of 
state of the inflaton during reheating in the minimal scenario. We illustrate this using chaotic 
inflation with a monomial potential as an example as in Eli. We then move on discuss how one 

^The upper bound may be lowered to 10® GeV in order to avoid gravitino overproduction if one considers 
supersymmetric models 0. 
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might use this theoretical uncertainty due to reheating to constrain the reheating temperature 
in some general universality classes of models in future precision experiments. In Section 
we then discuss a specihc example, the non-minimally coupled 0^ model, where the model 
predictions are robust and independent of the detail of reheating in the minimal scenario of 
reheating. Section is devoted to a short summary. 


2 Reheating corrections in single field infiation 

We consider a certain mode of the comoving wavenumber k which exits the horizon during 
inflation at /c = akHk- Considering the ratio of this scale to the size of the present horizon 
aoHo, we can relate it to various scales that correspond to critical moments in the cosmic history 
as [To] 


log 


k 

aoHo 


log 


I Qk ®rh Hk 
y fle ®rh Hq 


-Nk - iVi-h + log 




( 1 ) 


where the subscript e and rh denote, respectively, the end of inflation and the onset of radiation 
dominated era, viz. the moment of reheating. Nk = log(ae/afc) and Wh = log(arh/ae) are the 
e-folding numbers between the end of inflation and the horizon crossing of the mode of our 
interest, and the moment of reheating, respectively. Here Wh is highly dependent on the detail 
of reheating and is related to the reheating temperature T^h which appears implicitly in Orh/ao- 
Assuming the conservation of entropy after reheating, Orh/ao can be written as |H] 

flrh _ ^ ^0 


where gs* is the effective number of light species for entropy at the moment of reheating. 

Thus without detailed knowledge of the reheating process, any attempt to estimate Nk in 
Q is fundamentally restrictive even after we specify a model of inflation. In the following, 
we assume the minimal reheating scenario [S] where after inflation ends the inflaton coherently 
oscillates around the minimum of the potential, followed by decay to relativistic particles. If the 
decay is fast enough, Wh then depends only on how the inflaton oscillates around the minimum 
and is in principle also determined once a model of inflation is given. It can be written as [8] 


3(l + u;rh)^°®(prh) ■ 

Here, let us dehne = Prh/Pe < 1- This means that the reheating temperature can be written 
as 

T.h = -^PPe , (4) 

9* ) 

with g^ being the effective number of relativistic species at the moment of reheating, that is 
not necessarily the same as gs^. 
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Then, from Q Nk is estimated by 


Nk = - log 


k 


O^qHq 
1 — 3Wrh 

12(1 + Wrh) 


+ log ( ^ ) + log 


log (3. 



+ log 
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- 1/3 


TT 
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(5) 


For a fixed comoving scale k, the hrst and second terms on the right hand side can be very 
accurately determined by observations, whereas the third term and fourth term depend respec¬ 
tively upon the model of inflation and the underlying particle physics models after reheating 
ends. For realistic particle physics models, gg* ~ = 0(100). Given only depends on gg^ 

and logarithmically, it is accurate enough to take the general case where gg^, = g^ = 0 ( 100 ). 
Together they give the upper limit of Nk as long as <1/3. In general these contributions 
depend only weakly on the detail of the inflation model and the thermal history after reheating, 
as the dependence is logarithmic. The last term represents the correction from reheating. Once 
we hnd Nk, or (j)k precisely, we can accurately evaluate the spectral index and tensor-to-scalar 
ratio, including higher-order slow-roll corrections m 


nn-l = -6e + 27]+ [ 24a “ y ] - 2)e7/ -F + 


r = 16e 


1 + ( —4a H—— 1 c + ( —2a + x 1 h 



( 6 ) 

(7) 


where e = Tnl^iy'/Vf /2, g = mly" jV, = Tn^^V'lV^ and a = 2 - log 2 - 7 ^ 0.729637 

with 7 0.577216 being the Euler-Mascheroni constant. 

Let us hrst we consider simple power-law potential 

V(4>) = , ( 8 ) 


with a canonical kinetic term. Its effective equation of state during the oscillating phase can 
be expressed as [7] 

n — 2 

which is nearly constant if n is a global power index. The computation is straightforward even 
analytically PEI [ 12 ], and the model predictions of n-ji and r with the reheating correction are 
shown in Figure [Tj We can see that lower reheating temperature gives smaller uti and larger 
r, and that the uncertainty vanishes for V (0) ~ This is because for n = 4, from (§ we 
have tCrh ~ 1/3 and the last term in (|^ vanishes, thus the oscillation stage is indistinguishable 
from the radiation dominated epoch and the predictions of uti and r become robust. For an 
exhaustive study on other canonical single held slow-roll inhation models, see e.g. Pllsj. 

In fact, (|^ belongs to a universality class of models where the spectral index is inversely 
proportional to Nk, 

"K - 1 = -T . (10) 


where p = 0{1) is constant and depends on models, e.g. p = 2 for power-law inhation with 
n = 2 in (j^. Besides power-law inhation, this universal scaling relation also applies to a 
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Figure 1: riTi-r plot for the power-law chaotic inflation. The rightmost solid line denotes the 
instantaneous reheating immediately right after inflation ends at around Trh 10^® GeV. The 
leftmost dot-dashed line is the lower bound for Trh just above the nucleosynthesis bound 1 MeV. 
The intermediate dashed line symbols Trh = 10® GeV which is usually taken as the upper bound 
to avoid gravitino overproduction when supersymmetry is considered. Different power indices 
for the potential at n = 2/3, 1, 2 and 4 are marked by thick dashed lines. The Planck 2015 
contours (larger, TT-|-lowP) and (smaller,TT,TE,EE-|-lowP) are also shown. Notice that n = 4 
case gives a robust prediction, which, however, is excluded at 99.7% conhdence level by the 
Planck data. 


wide range of models HI including the inflation model [15]. With this universal scaling 
relation, we can ask what is the range we can constrain for W given a certain accuracy of 
measurements, which can be translated into that of the reheating temperature ATj-h- From 


(10), the uncertainty in uti can be related to that in 


AV, 


P 


An¬ 


il 


1 — 3rCrh 


{uTi-iy 12(1 -hWrh) 


log A/3, 


( 11 ) 


if we assume the main contribution to AN^ comes from the reheating correction term. With 
an accuracy of An-ji = O (10“^), which is expected from future observations such as the Euclid 
satellite 110. and a hducial value n-j^ = 0.968 from Planck[l], we see this gives A Afj. ^ 0.7p ~ 
1 — 2. Gonsidering a model that gives predictions within the observed range of utz, we deduce 
the uncertainty in Trh we can constrain is 


^ = o (10^ - lO^) 

Trh 


( 12 ) 


for canonical oscillation around a quadratic minimum, where Wrh ~ 0. For non-quadratic 
potentials the uncertainty can be different, and it can become as small as ARh/Rh = 0(10) 
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for Wrh ~ —1/3. On the other extreme case where w^h ~ 1/3, Nk is totally independent of T^-h 
and therefore Trh is not constrained from n-ji and r. These results are general as long as the 


model we consider are within the same universality class (10). 


We can consider as well a similar universal prediction on the tensor-to-scalar ratio, 


r = 


JL 


(13) 


where q and m are some constants of 0{1). For instance, quadratic power-law potential gives 
q = 8 and m = 1. The corresponding constraint on by the measurement on r is given by 


AW 


qljm /\ tp 


(14) 


For practical observability, we may consider only up to m = 2 that gives r = 0(0.01 — 
0.001), which is near the sensitivity of future observations on the S-mode polarization such as 
LiteBIRD |T7]. For m = 1, r is typically as large as 0(0.1). Thus with the prospect Ar = 10“^ 
we hnd AW = 0(0.1), e.g. for quadratic power-law potential AW ~ 0.45 [IH] so that we can 
constrain W with a good accuracy. 


ATrh 

^rh 


0 ( 1 ). 


(15) 


Even for smaller value of r with m = 2, AW/Rrh is better than or at least as good as what we 
can obtain from for r = 0.01, we can hnd ATi-h/W = O (10 — 10^). 

Similarly, we can do the same for the running of the spectral index, a^i. From (10), in 


general aji oc N, so that the constraint on W we can obtain from the measurement of is 


AW 


Aa 


7^ 


a. 


3/2 

7^ 


(16) 


with s = 0(1) being constant. Given the hducial value an = 10“^ and the accuracy ^an = 
O (10“^ — 10“®) that could be achieved in joint constraints using CMB and future large scale 
structure or 21cm observations like Square Kilometre Array |T9], we can expect AW/Rrh = 
0(10). Thus a detection of an in future experiments would give a better or compatible con¬ 
straint on Trh compared to that given by the measurements of nn- 

We should stress that in general the predictions of n:;^, r and a do not only depend on the 
inhation and reheating dynamics, but also on other phases after the end of inhation. Phases 
additional to the usual post-inhationary cosmic history where reheating is followed by radiation 
and matter dominated eras, if any, also induce uncertainties in W and thus can bias our 
predictions of the bound of the reheating temperature. We conhne ourselves within the minimal 
scenario and do not consider such scenarios here. 


3 A case study: non-minimal inflation 


The power-law potential we have discussed in the previous section, which belongs to the uni¬ 
versal class (10), is usually regarded as an ehective description at around the energy scale of 
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inflation. The effective eqnation of state ([^ is only valid when the power index n remains 
the same both in the large held region and the neighborhood of the minimnm. It is therefore 
interesting to stndy the extension to this monomial power-law inhation, like for instance the 
polynomial inhation stemmed from renormalizability [20] or some snperpotentials |2I]. These 


examples can be still described by ( 10 ) and after we simplify the polynomial to a piecewise 


potential, all the discussions are parallel. As in [ 20 ] , a natural cutoh of the power series will be 
at n = 4, which means the potential behaves as quartic in the large held region. It is already 
clear that the constraint on the reheating temperature becomes loose when n is approaching 
4, or, from another point of view, n k. A gives a robust prediction for 77 , 7 ^ and r independent of 
the detail of reheating. However, as now the large-held monomial quartic potential is excluded 
at 99.7% conhdence level by the Planck data [1], it is worthwhile to seek for other models that 
are also independent of the reheating temperature Trh. A natural consideration is to introduce 
a non-minimal coupling to gravity for the quartic monomial or polynomial up to quartic, which 
is also physically well motivated by the Higgs inhation model [22]- The action in the Jordan 
frame is 


S.i = / d^x 


m 


PI 


1 + e- 


m 


PI 


R 


\,rB,,pBA - 


(17) 


In Higgs inhation 0 is identihed as the standard model Higgs boson and there is a vac¬ 
uum expectation value v ~ 246 GeV in the potential corresponding to the mass rriH ~ 
125 GeV [23], which is neglected here. Performing the conformal transformation of the metric 
(1 -1- C0^/’^pi) Qijlvi we can work in the Einstein frame with the following action 


Se = / 


777 ^ 1 

~ ■j^s‘“B^xB„x - w'(x) 


( 18 ) 


where the canonically normalized held in the Einstein frame y is connected to 0 by 

^ ^ yi + (1 + 

d(j) l + 

and the ehective potential in the Einstein frame W is related to the Jordan frame potential by 


(19) 


W(4,) = 1 


4 (1 +fi^Vmpi)^ 


( 20 ) 


Note that since the conformal transformation factor to the Einstein frame is a function of 
0 only, perturbations and their correlation functions are equivalent in both frames [23] and 
calculations can be performed in the Einstein frame. For example, the power spectrum of the 
curvature perturbation is 


Rn 


A 1+ ^(1+ 0^ 


4 ’ 


7687r2 (1 + 

which can be used to replace A, and the energy density at the end of inhation as 

3 (^^1 +64^/3 + 128^-1)" 


( 21 ) 


de 

m- 


PI 8e (1 + 12^ + Vl + 64^/3 + 128e) ^ 


( 22 ) 
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In the minimal scenario of reheating (see also M), Wrh depends npon only the held dynamics 
near the minimum of the potential. For ^ 1, there is little deviation from the simple quartic 

power-law potential, so that tCrh ~ 1/3 and the predictions for n-ji and r are hxed irrespective 
of Ti-h- They he, however, on a curve rather than a single point on the n-ji-r plane depending 
on the value of On the other hand, for ^ ^ 1, there exists a critical value Xcr where the form 
of the potential changes m- 




for Xcr < X < Xe 


(23) 


-X 


for X < Xc 


where Xe ~ ''^Pi corresponds to the end of inhation and Xcr ~ \/273mpi/^ is the critical point 
where the dominant term in the potential changes. Thus, right after inflation the oscillation is 
sinusoidal as in a quadratic potential and w^h ~ 0. If the inflaton decays before the oscillating 
amplitude becomes smaller than Xcr, reheating happens during this effectively matter dominated 
stage. If the decay happens later, the potential is quartic and tCrh ~ 1/3. Again, it goes into 
the region where the prediction on n-ji and r is independent of the reheating temperature. 


We have to hnd the range of ^ where the approximation (23) is valid. We can see from (23) 


that if ^ is not large enough, the energy density at the critical value Xcr will be larger than 
that at the end of inflation and the potential exhibits only monimial quartic power, leaving no 
quadratic behavior as in (23). The value of ^ where this happens, can be found by solving 
p(Xcr(^*)) = Pe(^*) which, with and (1^, gives 


3^4 


1.457. 


(24) 


It is not much larger than 1 but we can still assume that the potential take the form of (23) 
when ^ 

Using the same procedure as in Section we can hnd W similar to (g. If ^ the 

predictions are identical to the instantaneous reheating, for which Nk is given by just (|^ with 
(3 = 1. If ^ the intermediate quadratic stage Xcr < X < Xe also affects the evolution and 
we have an extended version of (|^: 


Nk = - log 


k 


aoHo 

1 


+ log ( ^ 


+ log 


Hk 




1/4 


_3{1 + Wi) 3(l-Fta2)J \pej 12(1 +W 2 ) 


log ( ^ I + 


+ log 


1 — 3W2 


11 


- 1/3 


TT 

30' 


1 / 4 ' 


log/?, 


(25) 


where pcr is the energy density at the critical point Xcr, and Wi and W 2 denote respectively 
the effective equations of state for the oscillating phase before and after the critical point Xcr- 
Note that prh disappears in the hnal result as W 2 = 1/3, which means that all the predictions 
for a reheating temperature smaller than are the same. This is a strong assert implying 
the robustness of the predictions for the entire model as we can see in Figure The relation 
(25) can be extended to multiple stages, which can be used to study a piecewise potential for 
example including the hnal quadratic oscillation around the minimum at n ~ 246 GeV in the 
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Higgs inflation, which is also shown in Figure In general, a polynomial potential, if the 
characteristic scales for each term to dominate have a hierarchy structure as in (23), can always 
be approximated by a piecewise potential and studied similarly as in (25). 



Hr 


Figure 2: The predictions on the riT^-r plane of the non-minimal inflation model dl?] ). The 
rightmost dashed line corresponds to the instantaneous reheating with T^h 10GeV, while 
the solid line represents a huge range of the reheating temperature, Trh ^ min[^“^, 1] x 10^® GeV. 
This case is almost indistinguishable from the instantaneous reheating for ^ ^ 1, but approaches 
the predictions of the model for large If we consider a non-zero vacuum expectation value 
V, which we take v ~ 246 GeV, the predictions begin to change when the quadratic oscillation 
around v becomes effective, corresponding to T^h ~ 350 GeV. The leftmost dot-dashed line 
corresponds to the reheating temperature of 1 MeV, with a non-vanishing vacuum expectation 
value taken into account. The non-minimal coupling runs from 10“® (top) to 10® (bottom). All 
the predictions asymptotes those of the model when ^ is large, but the same point represents 
different e-folding numbers and different reheating temperatures in these two models. 

As we can see in Figure the predictions on n-ji-r plane are quite insensitive to T^h and 
are mainly dependent on the value of small ^ corresponds to the simple quartic power-law 
potential, while large ^ asymptotes the R^ inflation model. However, the same point in the 
asymptotic line predicted by the non-minimal coupling model with a large ^ and by the R^ 
inflation model corresponds to different reheating temperatures in these two different models. 
For example, the point in Figure]^ of {n'ji,r) = (0.9633,0.004024) represents a huge range of 
reheating temperature of 350 GeV < Tj-h ^ 10^® GeV in the non-minimal coupling model, while 
it only corresponds to a single but much smaller reheating temperature of 1 GeV in the R^ 
inflation model. It is becasue the effective equation of state is different in the two models. 
In the non-minimal coupling model, Wrh ~ 1/3 over a huge range of reheating temperature 
whereas rCrh ~ 0 in inflation. Thus although the non-minimal model provides the 1/N 
scaling to the tilt as R^ inflation in the regime we consider, the bounds on T^h are much worse 





here as we can see from (0- 

The uncertainty from reheating occupies only a very small region in the riT^-r plane when 
the potential is approximately quadratic. Note that the reheating temperature for the Higgs 
inflation lies in the region [2S] 

/ \ \V4 

3.4 X 10^3 GeV < < f — j x 1.1 x 10^^ GeV. (26) 


The largest possible uncertainty for e-foldings is given by 


AN, 


12 


log ( 


Pci 

\Pe 


< xlog^ ~ 2.28, 
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(27) 


for a typical value of ^ = 10®. Also, the corresponding n-ji and r with this uncertainty AN, 
taken into account are 


riTe = 0.9640 ±0.0007, , , 

28 

r = 0.0040 ±0.0003. ^ ^ 

An interesting point is that a lower bound of the tensor-to-scalar ratio r > 0.003714 is predicted 
irrespective of Note that a detailed study of reheating including the preheating stage in the 
Higgs inflation model is presented in |27] in parallel with |26j . 

4 Conclusions 

To conclude, we have studied the contributions of reheating to the predictions of n-ji and r. 
Assuming perturbative and fast reheating, the uncertainty from reheating is determined by the 
reheating temperature and the effective equation of state during the oscillating stage after 
inflation, which is determined by the inflation model. We have considered a universal class of 
predictions on the uncertainty in the reheating temperature, and found that future observations 
can pin down T^h within a factor of 0{1) optimistically if the oscillation around the minima 
of inflation potential is sinusoidal and Wrh ~ 0. We have also considered the inflation model 
with a scalar field non-minimally coupled to gravity, motivated by the Higgs inflation, as a 
typical example that has a different oscillation behavior. We have found that its predictions 
are insensitive to the reheating temperature, allowing precise predictions even without any 
knowledge of reheating. In studying this model, we have established an equation of N, for 
models involving different phases of oscillations with different equations of state, which can be 
easily extended to other cases like polynomial or piecewise potentials. 
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